Time-dependent mean field theory of tlie superfluid— insulator phase transition 
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We develop a time-dependent mean field approach, within the time-dependent variational prin- 
ciple, to describe the Superfluid-Insulator quantum phase transition. We construct the zero temper- 
ature phase diagram both of the Bose-Hubbard model (BHM), and of a spin-S Heisenberg model 
(SHM) with the XXZ anisotropy. The phase diagram of the BHM indicates a phase transition from 
a Mott insulator to a compressibile superfluid phase, and shows the expected lobe-like structure. 
The SHM phase diagram displays a quantum phase transition between a paramagnetic and a canted 
phases showing as well a lobe-like structure. We show how the BHM and Quantum Phase model 
(QPM) can be rigorously derived from the SHM. Based on such results, the phase boundaries of 
the SHM are mapped to the BHM ones, while the phase diagram of the QPM is related to that 
of the SHM. The QPM's phase diagram obtained through the application of our approach to the 
SHM, describes the known onset of the macroscopic phase coherence from the Coulomb blockade 
regime for increasing Josephson coupling constant. The BHM and the QPM phase diagrams are 
in good agreement with Quantum Monte Carlo results, and with the third order strong coupling 
perturbative expansion. 



PACS: 73.23.-b 73.50.-h 74.50.+r 75.45. -hj 
I. INTRODUCTION 

Phase transitions induced by thermal or quantum 
fluctuations have been studied in various mesoscopic 
systems. Examples are Josephson Junction Arrays 
(JJA) granular |^], and short-lenght superconduc- 
tors 1^ . Such systems have two different critical temper- 
atures Ti and Tq (Tq < Ti ) . Below Ti , they possess finite 
domains in which the electrons form the Cooper pairs: In 
each domain the condensate is described by the Cooper 
pairs' wave function ipj — A e**^ (A — beeing related 
to the pair density) and the system is globally resistive 
because of the absence of phase coherence between the 
Cooper pairs Below Tq, superconductivity may be ex- 
tended to a macroscopic scale (global superconductivity) 
if the system reachs the macroscopic phase coherence. 
If the characteristic energy scale of the system is much 
smaller than |A|, one can regard the Cooper pairs as true 
bosons and the global superconducting phase transi- 
tion can be studied by analyzing the critical behavior 
of strongly correlated bosonic models on a lattice [^,0. 
These exhibit two characteristic energy scales: the hop- 
ping amplitude t which accounts for the boson kinetic 
energy, and the Coulomb repulsion U which is the electro- 
static energy expense to make bosons spatially close. The 
BHM can describe the energetic competition between t 
and U . The global superconducting phase transition is 
controlled by the ratio t/U which is a measure of quantum 
effects. For t U, the strong quantum fluctuations of $j 
prevent the system from reaching the phase coherence for 
any value of the temperature. The condition t ^ U en- 



tails the classical regime: The system undergoes a phase 
transition at a finite value of temperarture Tq that be- 
longs to the ZJ-dimensional XY model's universality class. 
Below To the system is superfluid, while above Tq it be- 
comes resistive. If the Coulomb interaction is compara- 
ble with the kinetic energy then quantum fluctuations 
make vanish Tq, and drive the [Tq = 0) Superconductor- 
Insulator (SI) phase transition. The latter has been stud- 
ied in great detail both experimentally and theoret- 
ically [|6|. The superfluid phase is characterized by off- 
diagonal-long-range-order signalled by a non vanishing 



order parameter ^' 



(macroscopic quantum phase 



coherence). The insulating phase is incompressible and 
it is characterized by 4* = 0. In particular, due to the 
dimensional crossover, the SI phase transition belongs to 
{D -\- 1)-XY model's universality class for commensurate 
bosons' densities, whereas it is mean field-like away from 
such values . 

We recall that the number of bosons nj (at each site j) is 
standardly considered to be canonically conjugated with 
$ j . This establishes a competition between the quantum 
fluctuations of nj , and those of <I>j 's 

In a recent paper pO[ we formulated a Time- 
Dependent Mean Field Theory (TDMFT) of the BHM in 
order to investigate some aspects of the physical scenery 
just described. The TDMFT was based on factorizing 
slow/fast dynamics described by an effective form of the 
BHM Hamiltonian. The latter was derived within the 
Time-Dependent Variational Principle (TDVP) proce- 
dure, and relied on a picture of the system quantum state 
in terms of Glauber coherent states. In this approach the 
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hamiltonian degrees of freedom identified by construction 
with the parameters Zj (r) -r is the real time-of coherent 
states \zj), that is the expectation values {{zj\aj\zj)) of 
the operators aj {a~j) describing at each site j the de- 
struction (creation) of bosons |0 . 

We revealed that quantum effects concerning the compe- 
tition between the Coulomb term and the hopping term 
are embodied in the time dependence of the coherent 
state parameters Zj{T). The TDMFT involves a time- 
dependent, local order parameter which is assumed to 
represent the slowly varying part of Zj, and plays the 
same role of ijjj. In we have shown that ipj has a 
time-independent amplitude which is the analogue of A, 
and a time-dependent phase which is the analogue of . 
In particular, the phase's quantum fluctuations were de- 
scribed in terms of phase's time-fluctuations. 
The phase transition is signalled by a qualitative vari- 
ation of the time behavior of the local superconducting 
order parameter. In spite of the approximations involved 
by the TDMFT, indeed our phase diagram shows a good 
agreement with Quantum Monte Carlo (QMC) simula- 
tions Jl3t and Strong Coupling Perturbative Expansion 
(SCPE) Q. 

The purposes of this paper is to extend the TDMFT of 
the BHM developed in Ref. both to the SHM and to 
the QPM for constructing their zero temperature phase 
diagrams . To this end we establish a rigorous map- 
ping between the SHM, the BHM and the QPM based on 
the Holstein-Primakoff Realization and the Villain Real- 
ization of the spin algebra. In particular, we shall see that 
the existence of such a mapping is crucial to construct 
the phase diagram for the QPM within the TDMFT of 
the SHM. In this case in fact the explicit representation 
of the Hamiltonian in terms of coherent states is prob- 
lematic [0 due to the euclidean algebraic structure of 
QPM's degrees of freedom. 

In outline the paper is organized as follows. In Sec. 
II, we introduce the three models we deal with, we illus- 
trate some basic aspects of their algebraic structure and 
describe their qualitative phase diagrams. Based on the 
TDMFT, we construct in Sec. Ill the SHM's phase dia- 
gram. The latter will be shown to describe the SI transi- 
tion only in some interval of external magnetic fields. The 
discussion developed here will concern the SHM's phase 
diagram only for magnetic field in that range. The full 
description of the SHM's phase diagram will be reported 
elsewhere ITtI . In Sec. IV we employ the program devel- 
oped in Appendix A to recover the BHM phase diagram 
from the SHM one, and to obtain the quantum Josephson 
model's phase diagram from the SHM phase boundaries. 
In Sec. V we give our conclusions and further remarks. 
Appendix A is devoted to give a novel procedure follow- 
ing which both the BHM and the QPM are obtained from 
the SHM. In Appendix B we apply the TDVP method 
to work out the semiclassical BHM Hamiltonian and its 
dynamical equations. Then we formulate the TDMFT 



and employ it to construct the BHM's phase diagram. 
After reviewing the basic properties of spin generalized 
coherent states ||l^, the TDVP method is implemented 
in Appendix C for the SHM. In Appendix D we derive 
the phase dynamics of the QPM as a perturbation of the 
SHM minimum energy configurations. 

II. THE MODELS 

A convenient starting point for introducing models 
that exhibit the SI phase transition is the BHM [Q. 
It represents a boson gas of identical charges hopping 
through a Z?-dimensional lattice whose Hamiltonian 
reads: 

Hbh [^"^"j - 1) - ^] " ^ H + ' 

(1) 

where the operators := aja^ count the number of 
bosons at the site i, and the annihilation and creation op- 
erators fli, a] obey the canonical commutation relations 
[ui, aj] = Sij. The set {11, a| n^} is the basis generating 
at each site a Heisenberg-Weyl algebra /14. Also, the pa- 
rameters i, U of are the hopping amplitude and the 
strength of the onsite Coulomb repulsion, respectively, 
while the chemical potential /i fixes the average number 
of bosons in each site. 

The phase diagram of the BHM has been studied thor- 
oughly by means of mean field an d variational [|l^ ap- 
proaches as well as perturbative |14| and Quantum Monte 
Carlo techniques. At t/U — 0, the minimum energy 
configuration is characterized by an integer number n 
of bosons at each site, and a finite energy gap /i = 2U 
for the creation of particle-hole excitations. This reflects 
the Mott insulator (MI) behavior of such a phase which 
entails a vanishing compressibility. The MI regime sur- 
vives (except for the degeneration points with fi = 2nU) 
when t/U > 0, inside extended lobes attached to inter- 
vals /(n) — (2(n — 1)U , 2nU) of the fj,/U axis in the 
t/U-fi/U plane. Elsewhere, in the phase plane, the sys- 
tem exhibits a superfluid character, both compressible 
and independent from the filling. 

At the lobe boundary the appearence of the superfluidity 
is announced by the vanishing of the energy gap between 
the states corresponding to n (or n — 1) and n + 1 (or n) 
particles (or holes). Also, at the critical points the vari- 
ance I]^(<I') of the phase of the superconducting order 
parameter is reduced so much as the quantum coherence 
can take place. Indeed, the $i's quantum fluctuation 
survives also in the superfluid phase and they are rigor- 
ously vanishing (Sf($) = 0) only in the classical limit 
t/U — !■ 00. Such two features characterize the whole 
MI-SF phase boundary as well as the onset of superfluid 
state. 
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The QPM is deeply related to the BHM. It is largely 
employed for the description of quantum JJA's in which 
the phases $i of the superconducting order parameter 
are dynamically relevant, the fluctuations of the modu- 
lus A being negligible at low temperatures. Since among 
the islands of nanofabricated samples no Ohmic current 
flows and quasiparticle tunneling is actually negligible, 
the QPM Hamiltonian defined as 



QP 



(2) 



<i,j> 



with Qi := Ni — N^, can capture the physics of JJA. 
In (||), Ni counts the number of Cooper pairs which 
make the island deviate from the neutrality state with 
respect to the background charge. The standard as- 
sumption that phase operators $i are canonically con- 
jugated to Ni's, namely [$i,iVj] — 2ei Sij, requires co- 
herently that the Ni eigenvalues must range from — oo 
to -|-oo. The Coulomb interaction is described by the 
matrix Vij = 4e^(C~^)ij, where Cij is the inverse of 
the capacitance matrix. In the sequel we will assume 
V^j = Vdi^j with V = 4e^C^\ where Co is the self- 
capacitance. The external voltage Vx = N^/V enters via 
the induced charge eN^, and fixes the average charge on 
each island. The phase diagram of the QPM is similar to 
BHM's one (see j2^,|l). The MI lobes are attached to 
the interval I{Nq) = {{No - 1/2)V, {Nq + 1/2)V), where 
the average number of Cooper pairs is A^o = int{Nx). 
The degeneration points at E JV — are N^ — No/2. 
Outside the lobes the system is globally coherent and 
exhibits a superfluid character. 

The argument usually employed to map the BHM on 
QPM is heuristic: In the limit of large average num- 
ber of boson per site the bosonic field may be repre- 
sented as a,- ~ ,/n7e**^, a\ ~ e~**^yn7 (<i>j being her- 
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mitian), where the operators rij = Nj — N^ have neg- 
ligible fluctuations, and play the role of Cooper pairs' 
density. The QPM is recovered through the identifica- 
tions Nx = ^/ ^i.j ^n.'i = "-^ {'"^ the average 
boson density). Since the fluctuations of Uj are underes- 
timated in the QPM the superfluid region is smaller than 
the BHM one. 

It is worthwhile noting how a rigorous reading of this 
mapping leads to severe inconsistencies These are 

mainly due to the boundness from below of bosonic num- 
ber operators Uj that makes extremely complex the prob- 
lem of defining the hermitian phases canonically conju- 
gated to rij |2^. Nevertheless, the QPM can be defined 
in a consistent way since the operators Ni := —id^., and 
exp (±<I>i) with <I>i an hermitian phase are employed in 
formula (j^). 

The difficulties involved in mapping the BHM into the 
QPM in a direct way can be circumvented by using a 
third (spin) model as a bridge connecting the first ones. 
To this end we first describe the standard way to relate 



anisotropic spin-1/2 Heisenberg model with the physics 
of the SI transition, then we introduce the spin- 5' model 
that will be employed in the sequel. 
At very low temperature, few charge states are impor- 
tant if on-site Coulomb repulsions are very large. If the 
gate voltage is tuned close to a degeneracy, two charge 
states per island actually suffice to represent the rele- 
vant physics. Therefore the hard core BHM is equivalent 
to a spin-1/2 Heisenberg model (see Ref. Q) the model 
Hamiltonian of which is represented in terms of su{2) op- 
erators Sf, S^ , S~ in the fundamental representation. 
Substitutions Sf Ui — 1/2, S^ af , and Si — > al- 
low one to recover the hardcore BHM from the spin-1/2 
Heisenberg model. 

The zero temperature phase diagram of the spin-1/2 
Heisenberg model has been investigated in Ref. 24 
within the mean field approximation. This shows 
the presence of a phase transition between paramag- 
netic/Neel phases in which {Sf) t^O and {Sf) =0, and a 
canted state in which {Sf) = {S^) =0. 
Such magnetic orderings have a counterpart in the 
BHM's phases: the canted state (long range order in {Sf) 
and (S'f )) indicates superfluidity in the BHM; the param- 
agnetic/Neel phases (long range order in {Sf): \ {Sf) ^0) 
correspond to the Mi's (the role of fi is played by the 
Heisenberg model external magnetic field h). The repre- 
sentation in terms of the (anisotropic) spin-1/2 Heisen- 
berg model cannot work for the soft core BHM, this in- 
volving, in general, more than two charge states . 

We propose the spin-S XXZ Heisenberg model (de- 
noted so far by SHM) as a model capable of describing 
the physics of the SI phase transition when more than 
two charge states are involved. In this case the Hamilto- 
nian has the form of the anisotropic spin 1/2 Heisenberg 
model, but (the representation index) S > 1/2. It reads 



H,=J2iSt + S)[U{Sf + S-l)^h] 



E^ 



J2 (St Sr + S+ Sr) , (3) 



where we have taken into account only the on-site inter- 
action along z. In Appendix A we describe a procedure 
which maps the SHM onto the BHM and the QPM. 



III. TDMFT OF THE SHM 

In this section we construct the phase diagram repre- 
senting the SI phase transition for the SHM and compare 
it with that of the BHM represented in Figs. |l|,||. We 
develop a construction the main steps of which are very 
similar to the TDMFT of the BHM depicted in Appendix 
B; thus we will sketch them without comment when the 
analogy with the TDMFT of the BHM is evident. 
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The semiclassical SHM model is achieved by project- 
ing Hamihonian (^) on the su{2) coherent states (see 
Appendix B), thereby obtaining 



Semiclassicahy, E^(S'f) ^ 1 translates in L| 



Hs ^Y.^ULl - K)Ll - ^ + L*L,) + C, entails = V^^ _ ^2 

0j = Tp±/fi + (?!>j(0), then 



Within our scheme, the condition lA 
from the requantization procedure p6 
like variables L| (notice that 



= TO, Vj follows 
, p7t of the action- 
} = 5ij/K), and 
Since Eqs. (|l^) is solved by 



(4) 



where K ■= h - U{2S - 1), and C := Y.i S{K + C/(5 - 
1/2)]. The semiclassical equations of motion obtained 
within the TDVP method 121 read 



ihLj = \U{2L] + 25* - 



\)-h\L,+E,L]Y^U. (5) 

«e(j) 



Eqs. (^) are not integrable, since the constant of motion 
available are just two, namely the ^-component of the 
total angular momentum — 'Y^^ L^j and the energy it- 
self. The TDMFT decoupling simplifies Eq. (H) through 
the approximation 



L, L* 



M^L* + UM* 



MrM* 



(6) 

and entails 



which ensues from {Li — M.i){L* 
{Lj)r — M.j on long time scales (compare with (B7) 



and (BS)). The order parameter is defined as 



M \ 



(7) 



When (I) is inserted in (Q) together with the uniformity 
condition A4j = M. Vj, Hamiltonian (jj) reduces to a 
sum of on-site terms: TL^ Tij, where 

Hj = const + [U LI - h^)Ll 



(8) 



'i^{M,L*+M*L,~\M,\') 



The equations of motion derived from (S) simplify to 



± 



£/A-fp± 

qE.L^, ■ 



(11) 



We have parametrized the external magnetic field as 
A := h^^jU — 2m. The parameter (p_) is related 
to the choice {(f) j — aj) = {{4>j — ctj) = tt). The order 
parameter 



N. ^ ^ 



4m 
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has a vanishing (long) time average (Af) = 0, because of 
the time-dependent phase factor e*'^''±/''. Such a phase 
factor forbids the breaking of the original so(2) symme- 
try of the SHM in the xy plane. Thus, the oscillating 
behavior of Mj identifies the paramagnetic phase. No- 
tice that the condition (M) = can be realized also for 
^ 0: As in the case of the TDMFT of the BHM, 
the TDMFT of the SHM can describe the paramagnetic 
phase for i?^ > 0. 

The frequencies p±'s play the role of time correlation 
lenght governing the phase transition. The observations 
of the previous section as to the criticality of BHM can 
be extended to model (|3|): As expec ted, the critical ex- 
ponents z and V fulfill the same Eq. ( B16 ). The latter is 
left unchanged by the procedure mapping the SHM onto 
the BHM and QPM. Hence, this suggests that the same 
Eq. ( B16 ) governs the criticality of the QPM. 

Eqs (|ll|) play the role of the self-consistent equations 
of the TDMFT: They serve to eliminate the order param- 
eter from the energy (^) . The energy of the paramagnetic 
phase reads: 



ihLi 



(2UL' - h,) Lj ~ qE^L'Mj 



(9) 



and imply that Lz is no longer a constant of motion. In 
analogy to the BHM case of Appendix B, implement- 



ing the phase-locking condition 



Lj = |Lj|e* and Mj 



Uj — {0, tt} on 
a J e [0,27r]) 



successfully restores the basic feature dL^ /dt — 0. Due 
to Eqs. (^) the phases (pj obey the equation 



nLj4>j = (2[/L| - K) Lj - qE.L'^Mj 



(10) 



We examine first the spin dynamics related to the 
paramagnetic-like phase. Such a phase is identified by 
(Sj) sharply picked at one of its spectral values to: 
2^(5*^) <C 1; this condition (due to the uncertainty prin- 
ciple) induces strong quantum fluctuations of Sj and S'J 
which suppress the ferromagnetic order in the x—y plane. 



£,n{h,E^;p) = -U{A + m + S){m + S) 
2qEsm - p±-UA 



2qE^m? 



(12) 

[C/A + p±] {S'^-m^). 



The condition h^/U — 2m = (A = 0) identifies the de- 
generation points at which the canted phase reachs the 
axis E^/U = 0. The points on such axis correspond to 
static solutions of the equation of motions (^) which are 
obtained through p± = 0. Such a condition is the sim- 
plified form of the low frequency dynamics in the canted 
state. 

Now, we consider the fixed points of the dynamical 
equations (^) that identify the canted phase. In such a 
case the self-consistent equations read 



2UL', 



-L. 



(13) 
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The calculation of the energy minimimum of the canted 
phase is considerably simplified by applying the change 
of variable 



(14) 



which implies that the self-consistency equation can be 
written as 



(15) 



After noting that the order parameter Aij is a monotonic 
(increasing) function oipj we are allowed to eliminate the 
order parameter from the energy (|^) then minimizing it 
with respect to pj . The energy reachs its minimum value 
at Pj = 1 which corresponds to 



qE^ + 2U 



^0 , 



(16) 



where the index j has been dropped since the unifor- 
mity of the solution. We point out that inserting pj = 1 
in dTs] ) implies Mj — Lj on which the TDMF decoupling 
is based. The minimum energy is 



U 



= 2{m + S -1)- 



qE^m{'m + 1) 
U{2m+l) 



-1 



S ih + qSE, + Uf 



2{qE^ + 2U) 



(17) 



Energy ( |l7| ) gives the known (see |24|) classical value of 
the energy minimum (up to the constant US/2, see Ap- 
pendix B). Thus we conclude that the static solutions of 
Eqs. (^ correspond to the classical canted phase. 
Now we employ the energy (|l2|) to obtain the phase 
boundaries between the paramagnetic and the canted 
phases. The curves representing the m-phase bound- 
ary are identified by implementing the so(2) symmetry 
breaking through the limits p± — + and the vanishing 
of the energy gaps £„i — £m±i (compare with Eqs.(B21), 
(B22)). In particular, £m ~ £m+i provides the equation 



U 



2qE, 



= 



(18) 



where A_ = h/U — 2(m -I- 5) < 0, and := m(m -|- 
1)/(1 + 2m), Sra := [2(to -|- 5) + 1] (1 -h 2m)rl/ S'^ . A 
second equation is derived from £,„ = £m-i where A_ 
must be replaced by A+ = h/U - 2{m -I- 5) -I- 2 > 0. As 
in the bosonic case, the upper (lower) root issued from 
£rn — £m-i {£m+i ~ £m) givc the phasc bouudarics 



'1 - 



2U 



(2m + l)(2(m + S') + l) 



(19) 



IT 



2{m + S) 



qE^{m + l){m + 2) 
U{2m + 3) 



2U 



{2m + 3){2{m + S) + 3) 



(20) 



We study the phase diagram (g^, (|2C|) for m < -1/2. 
Such a restriction guarantees that the canted phase prop- 
agates up to Eg/U = at the degeneration points A = 
for any value of m < — 1/2 (otherwise the paramagnetic 
lobes meet one each other at A = 0, but at E^/U ^ 0). 
This feature is a common property of the phase diagrams 
which describes the SI transition. The above condition 
selects values of the external magnetic field h in the range 
—2U < h < {S — 3/2)U . Such a limitation will be relaxed 
in a forthcoming paper p^ ]. 

Fig. ^ shows the phase diagram of the SHM for S* = 10. 
We note that the coordinates {Eg/U)tr of the tricrit- 
ical points (the cusps of the lobes) are not monotonic 
as function of m (see also Fig. H) . In the lower part of 
the phase diagram they decrease monotonically (increas- 
ing m starting from its minimum value —S, h — —2 in 
Fig. ^ up to a critical value of m at which they are 
almost independent on m (we shall call "inversion re- 
gion" such a portion of the SHM phase diagram). In 
the lower part of the phase diagram, the shape of lobes 
is asymmetric. Instead, within the inversion region, the 
lobes are almost symmetric around odd-integer values of 



m + S*. Increasing S makes expecially wide the inversion 
regions. In such a region we can assume the effects of the 
L^'s dynamics independent on m. For larger values of m, 
{Eg/U)tr increases monotonically with m. Fig. |^ shows 
the non-monotonic behaviour of the tricritical points for 
various values of S. In line with the mapping developed 
in Sec. Ill, in the next section we show that the SHM 
phase diagram contains the BHM (QPM) phase diagram 
in the lower (inversion) region. 



IV. FROM THE SHM TO THE BHM AND QPM 
PHASE DIAGRAM 



In this section we apply the results of section HI to 
construct the phase boundaries of both the BHM and 
the QPM from the phase diagram of the SHM. 

We first recover the BHM phase diagram as a limit- 
ing case of the phase boundaries (p^, ( |2C| ) for 5 » 1, 
a ^ 0. The latter allows one to expand the energy of the 
paramagnetic phase ( p^ ) in power of a; upon relating the 
parameters E^, h with t, fx as in section HI A, we find 
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that the energy of the paramagnetic phase (|1J) reduces to the form 

P± + U{5 + 1) 



1 + 2 



qt 



[p± + UiS + l)]a + 0[a^] 



(21) 



Such a formula shows that, at the zero-th order in a , 
£m matches the on-site energy in the Mott phase ( B17 ). 
In Fig. H we compare the phase diagram of the BHM with 
the SHM phase diagram worked out using the energy jl^ ) 
for S = 55. 

Similarly, we can obtain the phase diagram of the QPM 
by taking the limit S", m ^ 1, <C 1 of the phase bound- 
aries (^9|), (pO|). As pointed out previously, we recall that 
a coherent states picture of QPM is not yet available in 
terms of coherent states of e(2) (which is the algebra of 
the QPM microscopic degrees of freedom, sec (A_3)) [ p^ . 
Hence, the direct application of the TDMFT, which re- 



2™t1 



US 



lies on the coherent states description of the hamiltonian 
operators, cannot be implemented in a direct way. The 
mapping outlined in section III B is crucial to bypass such 
difficulties. It allows, in fact, to construct the QPM phase 
diagram within the SHM lobe-like structure provided the 
values of m/S remain inside a suitable range (see be- 
low). As stated in Sec. Ill B, the condition S ^ |m| 
together with the identifications U = V, = /i*/2[/, 
and Ej = 5*^-5^(1 - /3'^) recast the SHM Hamiltonian in 
the QPM form (see appendix C). Correspondingly, the 
phase boundaries take the form 



UE, 



E^^T\ Ej+8 ^^2^(1-/3) 



(22) 



We note that for sufficiently large S small changes of 
m leave the parameter f3 almost unaltered. This is suf- 
ficient to make the SHM phase diagram periodic in Nx- 
The curves (|23) give a qualitatively correct QPM phase 



diagram for any P <^ 1. We fix /3 in ( ^21) to repro- 
duce quantitatively the QMC phase diagram. The Fig ^ 
shows the phase diagram of the SHM model for 5 = 55 
and — 0.028. The paramagnetic phases of the SHM 
turn in to the Coulomb blockade of the Cooper pairs. 
The canted state reveal the macroscopic phase coherence 
of the QPM. 



V. CONCLUSIONS 

In this paper we have been concerned with two different 
aspects of the SI transition. The former is the algebraic 
structure that characterizes the main models exhibiting 
the SI phase transition. The latter is the development of 
the Time-dependent mean field theory (TDMFT) for the 
spin 5* Heisenberg model with XXZ-anisotropy (SHM) 
that generalizes the approach previously elaborated for 
the Bose Hubbard Model (BHM). Remarkably, such a 
theory appears to be applicable to a large class of spin 
models. 

The analysis of appendix A concerning the algebraic 
framework in which the SHM is constructed, indicates 
that the correct way to map the SHM on the BHM is 
given by the Holstain-Primakoif realization of spin op- 
erators. Instead, the Quantum Phase Model (QPM) is 
related to SHM through Villain's realization of the spin 
operators. Thanks to the transparent geometric meaning 
of such realizations p5|,p7f both the BHM and the QPM 



can be issued from the SHM by considering the appro- 
priate sectors of the spin spectrum. To summarize: the 
BHM Hamiltonian emerges in the limiting case of spin 
vectors close to the south pole of the spin sphere, while 
the QPM behavior is found for spin vectors around the 
sphere equatorial plane. 

The TDMFT have been presented in an extensive way 
in appendix B, where it is used to investigate the BHM 
and its SI phase transition. Such a theory is based 
on combining the Time-dependent variational principle 
(TDVP), the coherent state picture of the model quan- 
tum dynamics, the Einstein requantization procedure, 
and a time-dependent generalization of the mean field 
decoupling. 

The central assumption in our theory is that at T = 0, 
time dependence of the semiclassical variables Zi{T) rep- 
resents the analogue of the quantum fluctuations of op- 
erators ai. 

Although Eqs. (B5) (the semiclassical counterpart of 



Heisenberg's equations for a^'s) have been simplified to 
describe an on-site dynamics concerning 2j's, they pro- 
vide a consistent description of the system's quantum 
phase transitions. This is due to the semiclassical re- 



quantization procedure (see (B14)), on the one hand, and 
on the time dependence characterizing the phase of the 
superconducting order parameter, on the other. The lat- 
ters make the order parameter phase the main respon- 
sible in driving the SI transition. Within the TDMFT, 
we have established a relation between the dynamical 
behavior of the (local) superconducting order parame- 
ter and the macroscopic phases exhibited by the BHM. 
The Mott phase has been shown to be characterized by 
time fluctuations of the phase of the superconducting or- 
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der parameter, whereas the superfluid phase is related 
to the static solutions of the (mean-field) equations of 
TDVP. The energy minimum coincides with the classical 
t/U ^ 1 superfluid case. 

Below, we compare the ordinary MFA (Mean-Field Ap- 
proximation) with the TDMFT. The first comment in or- 
der is that the time-independent MFA of the BHM [gj29| 
describes the SI phase transition by means of the sup- 
pression of the amplitude of the order parameter. The 
quantum fluctuations of the phase of the superconduct- 
ing order parameter play no role in the standard MFA. 
Moreover, the MI phase is characterized by an on-site 
energy where the hopping term does not contribute. 
The main difference between TDMFT and the standard 
MFA stands in the dynamical content of the definition 
of the superconducting order parameter: Within the 
TDMFT, the SF phase is suppressed by the (order pa- 
rameter's) phase^s time fluctuations just as phase quan- 
tum fluctuations destroy macroscopic quantum coher- 
ence. Also, in the MI phase only the time average of 
the local parameters ipj is vanishing (along large time 
scale). This does not imply that ipj is strictly zero inside 
the insulator as it happens in the standard MFA. 
The good agreement with Quantum Monte Carlo simu- 
lations and Strong Coupling Perturbative Expansion is 
a posteriori confirmation that the superfluid phase is al- 
most classical: Quantum fluctuations are strong in the 
insulating phase; they drive the SI phase transition, and 
are negligible in the superfluid phase, except in the re- 
gions very near to the phase boundaries. 

Based on the introduction of the spin coherent states 
performed in Appendix C, the TDVP effective picture of 
quantum dynamics and the TDMFT of Appendix B have 
been extended to the SHM in Sec. III. Our analysis has 
revealed a quantum phase transition between a param- 
agnetic and a canted phase. The phase diagram exhibits 
a lobe-like structure. Inside the lobes the phase is para- 
magnetic; elsewhere the system is in a canted state. The 
energy in the canted phase (represented by the station- 
ary solutions of the equations of motions) coincides with 
the known classical energy. 

We point out how the ordinary MFA (linearization of the 
xy exchange term in (^) leads to study the reduced (on- 
site) Hamiltonian in the su{2) enveloping algebra (due 
to the (S"^)^ term) which prevents one from diagonal- 
izing the Hamiltonian for generic S |Q. On the other 
hand, fixing the representation index leads to study the 
matrix-Hamiltonian whose rank increases with S. Thus, 
studying the spectrum of (|^) for sufficently high S (which 
is done in the present paper) would be very problematic. 
Otherwise the classical analysis (see f.i. [^) cannot cap- 
ture the lobe-like structure of the SHM phase diagram. 

The SHM's phase diagram contains the BHM's and 
QPM's ones. Using the strategy developed in Sec. Ill we 
have recovered the BHM energy as well as its phase di- 
agram. As we pointed out previously, implementing the 



TDMFT to the QPM is problematic since the semiclassi- 
cal description of the model in terms of the e(2) coherent 
states is not available. Nevertheless, the phase diagram 
of the QPM has been obtained as an appropriate limit 
(see appendix A and IV) of the SHM phase diagram (the 
lobe structure is (locally) periodic around m for changes 
TO -I- Am such that ISm/S <^ m/S). This picture ap- 
pears to be consistent with the QPM form assumed by 
the SHM (effective) Hamiltonian close to its ground-state 
configuration (see Appendix D), and suggests that pos- 
sibly the (purely quantum) SHM can be reduced to the 
QPM form as well. 

Other perspectives are opened by the present study. 
First, the tools we used to map the SHM on the BHM 
and the QPM show how the enveloping algebra charac- 
terizing the BHM and QPM have common root in the 
enveloping algebra of the SHM. The mapping between 
these models can be seen as contraction procedure of 
the enveloping algebras underlying the BHM, SHM and 
QPM [Hj. This suggests the fact that the universality 
class might be preserved by contraction. 
Second, since the equations of motion (of the BHM and 
the SHM ) obtained by the TDVP have been considered 
within the simplified form entailed by TDMFT, the dy- 
namical approach refined out the TDMFT should imply 
a more accurate description of the superconducting order 
parameter dynamics as well as a better understanding of 
the low temperature regime. As to this point, other im- 
provements can be achieved by constructing trial wave 
functions able to account more accurately the micro- 
scopic physical processes (see [^). Finally, the TDMFT 
succeeds in capturing the main features of the quantum 
phase diagram of a rather large class of models. Promis- 
ing applications of the TDMFT to other systems are ex- 
pected due to its feasibile character. 
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APPENDIX A: MODEL MAPPING 

In this appendix we give a new proof of the equiva- 
lence between the SHM, the BHM, and the QPM. The 
key operation of such a mapping consists in studying the 
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model (H) for sufficiently high values of representation 
index S, and in using the Holstein-Primakoff Realization 
(HPR) and the Villain Realization (VR) of the spin al- 
gebra. 



is consistently reproduced at the classical level. We shall 
see in Sec. VI that the SHM phase diagram matches the 
BHM phase diagram for a ^ 0. 



1. From the SHM to the BHM 

We rewrite spin operators in Hamiltonian (H) by means 
of the HPR of the spin algebra su{2) 



S+ = V2Sal^l~n^/{2S) 
Sf = n,- - S 



(Al) 



and note how, in view of the formulas of Appendix B, 
the eigenvalues of the secular equation for the i-ih spin 
{SI + 5)|m„ S) = (to, + 5)|to„ S) ranging in {0, ... , 25} 
identifies with the eigenvalues of n^. Consistently, the 
spin states \i!ni,S) coincide with the number operator 
eigenstates jn^) up to the reparametrization Ui = rrii + S. 

The condition m/S = {rrn + S)/S ~ a allows one 
to obtain Hamiltonian (|l|) from (||). In particular, the 
xy-exchange term of gets in the hopping term of (|l|) 
with the hopping amplitude t ~ 2(1 — a)SE^ — + 2SEg 
for S —^ oo, and rii finite. In the same limit, the spec- 
trum of + S = rii ranges in {0, . . . , oo} thus repro- 
ducing the spectrum of bosonic operators n.i. The z- 
antiferromagnetic term and the coupling with the exter- 
nal magnetic field term becomes the Coulomb interaction 
and chemical potential terms repectively {h = fi). 
Also, since a = {mi -\- S)/S {a has to be viewed as an 
order of magnitude independent on the site label i), the 
above limit corresponds quantum mechanically to select 
spin states close to lowest weight vector jS", — 5) of the al- 
gebra sui{2) Q at the z-th site. In this respect, we recall 
that, within the HPR, jS", -5) = |0), where a^lO) = 0. 
The effect of taking a ^ 1 is illustrated by means of 
the semiclassical spin vector L (defined in Appendix B) 
lying on a sphere of radius S. Such a limit leads to se- 
lect those vectors neighbouring to the south pole of the 
sphere. Consistently, spin coherent states (see defi- 
nition (pl[)), that pertain to the Hilbert space of bosons 
within the HPR, can be shown to tend to the bosonic ones 
\zi) for a ^ (see (B2)), while the same limit makes the 
spin vectors on the sphere coincide, via stereographic pro- 
jection, with the points of the south pole tangent plane 
iQ. These are in bijective correspondence with the com- 
plex numbers Zj labeling the bosonic coherent states. In 
fact, formulas (C3), ( |C^ ) clearly show that L| — 5 is 
achieved for negligible j^^ p which also entails L* 2S£,j, 



L 







2Sq. 



This in turn leads to the identifications 



Lj/V2S = V2Sa, and 



z* = L*/V2S = V2S^j,z,i 
makes coincide brackets (pT]) with brackets (B6). Hence 



2. From the SHM to the QPM 

We write the spin operators in Hamiltonian ^ by 
means of the Villain realization of su{2) This is 

based on the formulas 



^^(5 + 1/2)2 -(5| + 1/2)2, 



(A2) 



Such operators fulfil the su{2) commutation rules pro- 
vided the action-angle operators Sj, satisfy the e(2) 
commutators 



[Si 



-,±i4 



±6j,ie 



±i4 



. (A3) 



The QPM (|) is obtained as a limiting model from 
Hamiltonian (^) when considering the first order in 
(3 = \m,i\/S ^ 1 for S ^ 1. The ferromagnetic part 
of (0) reduces to the Josephson term with coupling 
E, = {S + 1/2)2^:^(1 - /32) which in the limit S ^ oo 
becomes Ej = S'^E^. As in the case of a, the parame- 
ter (3 = \mi\/S must be regarded as a site- independent 
order of magnitude. Consistently, the spectrum of Sf 
will range from — cx3 to +oo thus reproducing the spec- 
trum of the unbounded operator Ni in (||). The rest of 
Hamiltonian (^) maps to the charging term of (|^) pro- 
vided V =U, and N^ = [h + U{1 - 2S)]/2U. It results 
Hs -^QP + C'o, where Co ^ -{U + hf/W. 
Geometrically, (3 amounts to select vectors neighbor- 
ing the equatorial xy-plane of the semiclassical sphere of 



the limit described above for the hamiltonian operators 



radius S. This is well illustrated via formulas (C2), (C4) 
that provide |6P — 1 ^ 1 (l-^fl S) as counterpart 
of the above condition (3 <^ 1. Consistently, the sphere 
equation {L^f + \L*\^ = 5^ shows that \L*\^ ~ S"^ . Since 
both the Josephson coupling and the hopping amplitude 
contain the factor E^, the formula E, =tS/2 holds for 
1<5,Q!<1,/3<1. Thus S plays the role of the boson 
density n (see Sec. II). We point out that the effective 
hopping coupling in the BHM is reduced by a factor S 
when compared with the Josephson coupling in QPM; 
consistently, the SF region of the BHM's phase diagram 
is smaller than the SF region of the QPM's one. 

We note how, when considering the perturbative dy- 
namics around the ground-state configuration of the 
SHM (see Appendix C), one obtains a QPM- like behavior 
(that is, having a pure phase's dynamics) without per- 
forming the limit S ^ oo. We shall see that the phase 
diagram of the QPM is obtained for /J^ <^ 1 [Eel . 
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APPENDIX B: THE TDMFT OF THE BHM 

In this appendix we apply the TDVP (see Refs. ||2^ ] 
and fl^ ] for a general review) to the quantum dynamics 
of Hamiltonian (|l|) , and implement the TDMFT for the 
BHM. 



1. Time-dependent Mean field theory 

The initial step of the TDVP method amounts to 
finding a solution of the Schrodinger problem {ihdr — 
H)\^) = by approximating the exact (unknown) 
solution 1 4") through a macroscopic state |$) whose 
time evolution is imposed to obey the weaker form of 
Schrodinger's equation {^KiUdr — ^^)|*i>) — 0. Upon set- 
ting |<I>) = exp{iS /fi)\Z) one obtains 



S = ih{Z\dr\Z) ~HiZ) 



(Bl) 



(Ti-iZ) = {Z\H\Z)) which represents the key equation of 
the approach. 

The building blocks of Hamiltonian (|l|) are operators of 
the A^g-site Heisenberg-Weyl algebra hi{Na) — {I, a.i, a\ : 
i E A}, Ns is the number of sites of the lattice A, but 
actually belongs to the enveloping algebra A of W{Ns) 
because of the quadratic terms n^. This motivates the 
choice \Z) :— ®j\zj) as the trial macroscopic state, en- 
tailing 



(B2) 



where the states \zi) are the Glauber coherent states ful- 
filling the secular equation ai\zi) = Zi\zi) for the boson 
lowering operator U j, a t each site i. In this case the ef- 
fective Lagrangian (Bl) becomes 



S[Z] =i'h'^ '^{ziZi - ZiZi) - Ti{Z) , 



(B3) 



where Ti-iZ) = {Z\H\Z) -the semiclassical model 
Hamiltonian- is given by 



y)\zi 



2 {ziZj + ZjZii 



(B4) 



The equations obtained variationally from ([B3|) 



ihzi = ^iizi + 2Uzi\zi 



IE 



(B5) 



J6(i) 



account for the dynamics of variables (expectation val- 
ues) Zi = {zi\ai\zi). Eqs. (B5) describe an hamiltonian 
fiow in that they can be equivalently obtained through 
the standard formulas ihzj = {z^, 7i}, where the Poisson 
Brackets 



{fiZ,Z),g{Z, Z)} 



-y 



d£dg_ _ dg df 

dzj dzj dzj dzj 



(B6) 



specifically, {zk,Zj} — Skj/ifi have replaced the basic 
commutators [a.i,aj] — Sjj within the TDVP semiclas- 
sical framework. Eqs. ( |B^ ) are not integrable, since the 
only known constant of motion, apart from the energy, 
is Af = J2i l-^iP; the semiclassical version of the total 
particle number N = "^i^i. The presence of the non- 
linear {/-dependent term prevents one from decoupling 
them in the dual lattice space. 

The TDMFT procedure is, in a sense, the analogue in 
a dynamical contest of the Mean Field Approximation 
(MFA) usually employed in Statistical Mechanics and is 
based on a well known microscopic picture of superfiuids 
illustrated, e.g., in Ref. |^^. It leads to simplify the struc- 
ture of Eqs. (|B5| ). We set at each site Zi — ^pi + rji, where 
tpi is a slow variable, whereas rji is a fast oscillating term 
which describes the high-frequency part of the dynamics 
taking place on the hopping interaction time scale. Also, 
we assume that (z^ — ipi){zj — — ViVj ~ 0- Thus 



of the bj's 



when the time scale r is larger than that 
{{•)t denotes the time average). Such time 
averages coincide with statistical averages (in the Gibbs 
ensemble) under the ergodic assumption. The onset to 
the (macroscopically) ordered phase reflects the presence 
in the system of stable, slowly varying components of the 
lattice dynamics corresponding to the ipj^s. This means 
that any Zj is strongly attracted to its average value ^pj 
(namely that the collection of ijjj 's defines the dynamical 
system's attractor). Dynamical regimes where the long 
scale-time behavior of Zj is not described by an asymph- 
totic slowly varying function ^pj is related with the dis- 
ordered phases of the system. The above considerations 
imply the TDMF decoupling 

ZtZj= {zi - 1pi){Zj - Ipj) + ipiZj + ijjjZi - Ip^lpj 

« Tp^Zj + {pjZi - ip.ijjj . (B7) 

The dynamical scenery just depicted leads thus naturally 
to defining 



(B8) 



as the macroscopic order paramenter revealing when or- 
der issues from the the lattice dynamics. Using formula 
in a. modifies the kinetic term as follows 



ZiZ'\ -\- Zq Zi 



(B9) 



where ipj = ipi for j G {i) (smoothing condition). The 
resulting Hamiltonian reduces to the decoupled form 
T-Lmf = J^jT^j^ where 
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= u\z, 



qt 



(BIO) 



and exhibits a dimensionality dependence entering via 
the numbers of nearest neighbours q. The (decoupled) 
equations of motion ensuing from Timf read 



(Bll) 



When compared with the exact ones, (B5), they im 



ply the relation qipi « X]je(i) -^j consistently leading to 
an identity once the time average is carried on and the 
smoothing conditions is used. A further effect coming 
from the linearization ( p39| ) consists in the fact that the 
total particle number Af = J2i does not have any 
longer vanishing Poisson brackets with Timf- Restoring 
such a basic feature is performed by considering Zj with 
an appropriate time dependence. To this purpose we look 
for solutions of Eqs. ( [Bll| ) where 6j, Xj, the phases of 



^3 = IHe'""' 



(B12) 



respectively, obey the phase locking condition {9j — 
Xj) = {0,7r}. Then M is constant due to the fact that 
d\zj\'^/dt = iqt{zj'ifjj — Zjil!j)/2 = 0. Moreover, the fur- 
ther condition d\ipj\/dt = consistently makes Ti-mf con- 
stant as exp ected for the total energy. 
Due to Eq. ( Bll ), the phase 9 obeys the equation 



h\z,\9, ^ {2U\z,\' - ti)\z,\ 



(B13) 



where s = ±, depending on how the phase locking con- 
straint is implemented. In spite of its simplicity, such an 
equation is able to characterize both the MI phase and 
the superfluid phase in terms of phase dynamics. 

We examine first the dynamics related to the MI. In 
this case, ipj must have a zero time average along macro- 
scopic time scales. Such a behavior occurs when the uni- 



form filling conditions nj 



for all i (we identify here 



number operators rij's with their integer spectral values) 
is inserted in (B13) by setting 



= n e TV. 



(B14) 



Such a substitution is the natural consequence of the re- 
quantization process ||2^,|2^ of the action-like variables 
|2jP (notice that {|2.i| ,0^} — dij/fi) strongly requested 
from the pure quantum character of the MI. At t/U ~ 0, 
where the system is integrable (since it reduces to a set 
of uncoupled, nonlinear oscillators) indeed 9j and jz^p 
represent the pairs of action-angle variables of the sys- 
tem. For small values oi t/U such a feature is still true 
as a consequence of the fact that the nonlinear oscillators 
are weakly interacting in the MI regime. Hence, in the 
spirit of Einstein's requantization procedure (see [p6[). 



their orbits are still homotopic to those of the integrable 
case which entails again \zj 



Eq. ( B13D is easily showed to be solved by 6'j(r) 
X±T/h + 9j{0), with X± defined through 



|^,|=±^(A±-C/5) 



(B15) 



Here we have parametrized the chemical potential as 
S = ^/U — 2n, and A_ (A+) is related to the choice 
9j — Xj = TT i9j — Xj — 0) (notice that the index 



j does not label X± since the request 



leads to — ^Jn at each site). In the present the- 
ory, the frequencies A± play the role of time correlation 
lenght governing the phase transition. Our theory gives 
A± = U^/n{^ — ^c) for fixed t and A± — q\^i\/2{t — tc) 
for fixed fi (fic and tc are the critical values of /i and 
t). Upon defining the critical exponents z and v as in 
Ref. 0, we argue that g 



zv — 1 



(B16) 



By replacing in the reduced Hamiltonian ( B1C| ) the value 
of IV'il provided by Eq. (BK ), the energy of the MI reads 



EylpL,t;\±)=n 



qt 



{\± + U5f + U{5-n) -2A± 



(B17) 



where the subscript n reminds us that the filling n 
is accounted for. The oscillating behavior of = 



iOj (0) 



having a vanishing long time av- 



erage, identifies the MI. This, in fact, implies that the 
gauge symmetry breaking expected in the SF phase can- 
not take place. Notice that the ordinary (time indepen- 
dent) MFA cannot describe the MI for i > 0, since the 
hopping term of the reduced Hamiltonian is canceled by 
the vanishing of the order parameter, ip — Q. Within 
our scheme, instead, the condition {^)r = can be re- 
alized also for -0 7^ 0. The degeneration points selected 
by = 2n (i.e. 5 = in I{n)) are extreme limit- 

ing points for which the superfluid phase is extended up 
to t/U — 0. They will be identified with the meeting 
points of the lobe boundaries. Such points characterize a 
static phase due to A± = (see Eq. ( B13 ). We interpret 
the stationarity which distinguishes the solutions of the 
semiclassical equation of motion as the trait characteriz- 
ing the SF phase in which Y?{e) (classical SF). This 
is but the oversimplified version of the low frequency dy- 
namics expected in the SF phase that should correspond 
to the condition S^(6') < 1. 

Let us consider the fixed points of dynamical equa- 



tions (Bll) that, as we concluded above, identify the clas- 



sical SF phase. Such solutions (the trivial case ij = due 



to Zj = ■tpj = is excluded) allow us to recast Eq.(Bll) 
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in the form ipj =2[{2U\zj\'^ —^)/qt\zj making ipj a func- 



tion of Zj. Inserting ipjlzj) in (BIO) reduces the energy 



associated to the Hamihonian Hi to 



:(^,t,Zj) = |Zjf 



-{^i-2U\z,ff+^l-3U\z,\' 



(B18) 



The quantity e{^,t,Zj) is the on-site energy accounting 
for the absence of dynamics. The hmit A± — s- 0, in fact, 
shows that En{n,t; X±) e{^,t,Zj) provided n = 
The lowest value of energy ( p318 ) and the value in- 
volved for Zj are obtained by minimization. They are 
given by 



-U\z,\^ 



AU 



(B19) 



tips are singular points of the energy. Their coordinates, 
obtained imposing fi-{t) = fJ'+(t), reads 



tc = U /qn 



(B24) 



and ^{tc)/U = 2n — 1 — (l/2n). The values of tc ob- 
tained within the present theory, can be compared with 
QMC and SCPE |4| (Figs. |l|,|). Contrary to the result 
obtained in Ref. [H, our phase diagram has a concav- 
ity independent of the dimensionality. In ID we find a 
good agreement with QMC and SCPE. Upon recalling 
that our construction relied on Eq. (Bll) -this incorpo- 
rates the time-dependent mean field approximation- it is 
important to note that the concavity of lobes in D > 2 
might be improved by implementing the requantization 
process directly on the Eqs. (B5). 



\z4' 



(m + g^/2) 
2U 



(B20) 2. Remarks on the SF ground-state 



respectively. The phase of z» can be set to zero since 
the gauge symmetry breaking characterizing the ground- 
state configuration. It is worth noting that inserting 
in the expression 'ipj(zj) implies that ipj = Zj so that 
the minimum energy configuration naturally fulfils the 
consistency condition on which our TDM FT is based. 

Now, we employ the expression (B17) for the on-site 
energy to construct the BHM's phase diagram. In the SF 
phase, the states with n and n -I- 1 (adding a particle) , as 
well as the states with n—1 and n (adding a hole) must be 
degenerate. The curves representing the n-lobe bound- 
ary are identified by implementing both gauge symmetry 
breaking through the limits A± — > and the vanishing of 
the energy gaps En — En±i- In other words we require 



hm {En 



En^ 



= 



lim (En-i - E,,) = 



(<5<0) 



(<5>0) 



(B21) 



(B22) 



For solving Eqs. (B21)-(B22) we introduce the variables 
5± = IJ-IU - 2n + (1 ± 1). By inserting (5+ > {5_ < 0) 
m eq. ( p2ll ) ((p2^)), and defining r = qt/4U, one gets 
the quadratic equations (5^. -|- 2r5± — {2n =F 1) = 0, that 
furnish the pair of two-branched curves 



U 



2n-l =F 1 



It 



(2nTl) 



As to the effectiveness of the approach just illustrated, 
two important observations are in order. The first is 
that our finding ( |B19|) , ( B20 ) concerning the ground- 
state configuration is remarkably confirmed by two other 
procedures. The other concerns the macroscopic phase 
at the phase transition (see Eq.(B27),(B28)). 

The minimum energy of the SF phase can be calcu- 
lated from (B4). It represents the exact value of the BHM 
ground-state energy in the classical hmit t/U » 1. The 
minimum energy is readily obtained by rewriting first the 
hopping terms as {ziZj+ZjZi) = \zi\'^ -\-\zj\^ — \zj—Zi\'^ , and 
by noticing then that the choice = ^ for each site en- 
tails the lowest value of the hopping term since the only 
positive contribution \zj — Zi\'^ vanishes on each bond. 
Upon minimizing the resulting expression of the energy 



7i* 



'-lie 



(B25) 



by setting = 0, one obtains = {fi+qt/2)/2U 

that matches exactly (B19): The ground-state energy 



(B18) coincides with the minimum of the exact (semi- 



classical) energy. Our approximation scheme thus repro- 
duces the correct value of ^ as well as the corresponding 
value of the ground-state energy. 

The ground-state energy eigenvalue can also be obtained 
once the Hamiltonian operator (^ is linearized via the 



(B23) standard procedure nj « 2w 



which is realiable 



for t ^ U. This yields the diagonalized Hamiltonian 



The lower branch n+(t) and the upper one n~{t) consti- 
tute the boundary encircling the n-th lobe. The substitu- 
tion of /i-(t) and /i+(i) in (B21) and (B22) respectively, 
provides the on-site energy values involved in the two 
cases, namely En-i{fii,t) = i?„(/ii,i) = Un{n — 1) and 
Enip-2, t) = En+i{fJ.2^ t) = Un{n -\- 1). The two branches 
are therefore separated by an energy gap. Thus, the lobe 



~ -[/iV,i/2 +^[2i7i^ - /i - tg{k)\ h+bk , (B26) 

k 

when the operators hk = Ns^^^'^TijUj eyi^[ikj] of the k 
modes are used. We have introduced g{k) — T,rCos{kr) 
[r G [^,D] on a D-square lattice) where kr is the r-th 



component of k. Hamiltonian (B26) clearly shows that 
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its lowest eigenvalue is obtained through the depletion of 
any mode fc 7^ (boson condensation in the state with 
fc = 0). As a consequence of the consistency condition 
(n) = V (=> N = Nsv), once more the energy is mini- 
mized hy v={fi + qt/2)/{2U). 

A comparison with the quantum ground-state energy 
-known exactly in the case t/U = 0- is important as 
well. For n/U & I{n) the eigenvalues of H^^ with inte- 
ger filling n = N/N„ Ei{nj}) = T,j[Un^. - {p + U)nj] 



reach their minimun value i?* = —UN^ 



for 



and ^/U = 2n—l. The on-site energy e* ~ — (see 



( B20D ) is found to attain exactly its quantum counter- 
part E^/Ns = -Uri^ in the limit t/U 0, fi/U 2n, 
namely at the point of I{n) representing its top. 

The second observation pertains the action S that rep- 
resents the phase of the macroscopic state | $) . It raises 
a special interest since it is itself a macroscopic quantity 
and thus is viable to experimental observations. In the 
following we compare S with Smf (i.e. S in our TDMFT) 
as well as Smf in the Ml with Smf in the SF. 
When Eqs. (B5) are inserted in (^^ ) then S reduces to 



S = U^,\z,\ 



(B27) 



where the explicit form of Zj(T) is known only once 



Eqs. (B5) have been really solved. 



Inserting Eqs. ( Bll ) instead of (B5) in (BS), and replac- 
ing Ti. with its mean field version Ti-mf , involves 



s 



mf 



[U\z,\^+t{z,^,+^,z,)-2t\^,\' 



(B28) 



If Zj « Vii then Smf and S have essentially the same 
form. Such macroscopic quantities may actually co- 
incide at the low temperature regime if the dynamics 
of both Ti-mf and Ti. have solutions characterized by 
\zj{T)\'^ « const. Within the present TDMFT Smf = S 
since the condition rjifij w implies Zj « ipj; further- 
more, |2;j(T)| is strictly time-independent. 

In the Ml the requantization rule ( B14) m ust be used. 
The frequency Smf is obtained from (828) by insert- 
ing ( B14 ). Smf reads as 



Smf = UNs' 



(B29) 



A transition that changes the filling from n to (n + 1) 
involves a change of the phase frequency amounting to 
UNs{2n -1-1). The action density in the superfluid phase 
reads: Smf = UNs\z\'^, which compared with the corre- 
sponding formula ( B29| ) shows that the frequency is not 
quantized. In the SF phase transitions between differ- 
ent configurations (different values of the filling) occur 
continously. 



10 ■.= Ds{i)e^^^\S,-S), 



(CI) 



where is the raising operator of the angular momen- 
tum algebra su{2) generated by 6*2, Sx = (5+ + S~)/2, 
and Sy = {S^ — S~)/2i which fulfil the standard com- 
mutators 



[S',S^]=±S^2 , [S+,S- 



2S, 



(C2) 



Ds{(,) = 1/(1 + ICP)'^ represents a normalizing factor, 
whereas jiS, — S"), the so-called maximum weight vector, 
satisfies the equation S'~|S', —S) = 0. The action of 



S^\S,m) = y/{STm){S±m+l)\S,m±l) 

is represented on the standard basis {\S, m); \m\ < S} the 
vectors of which obey the secular equation S^\S,m) = 
mis', to). Making explicit the action of 5"+ in (CI) sup- 
plies the spanned form 

s 

lO^DsiO E Cm{S)r-''\S,m) 

m=-S 



of 10, with Cm{S) = ^/i2S)\/{S -m)l{S + my.. Based 
on such formula one is able to calculate the expectation 
values 



s 



1 



L* = {S+) = S 



2e 



1 



(C3) 



(C4) 



where (•) = • and therefore to reconstruct the 
sphere equation (L^)^ + {L^)"^ + {L^Y = S^: where 
L = L"" - iiy = (L*)* for the classic spin {L'',Ly,L^). 
Such an equation, in turn, can be viewed as the classic 
counterpart of the (quantum) Casimir equation (S"^)^ -I- 
(5^)2 + (5^)2 S{S+1) for (S'^,S'S',5'") in terms of SCS 
picture and illustrate the semiclassical content thereof. 
Upon introducing the macroscopic wavefunction 



1$) ^ e'^/^\C) 



(C5) 



where the trial state - = ®i and |^i) is the SCS for 
the i-th spin {Sf , Sf , Sf ), then one easily constructs the 
TDVP semiclassical dynamics relative to spin Hamilto- 
nian (H) 



APPENDIX C: SPIN COHERENT STATE 
PICTURE OF XXZ MODEL 

Spin coherent states (SPS) \£_) are defined as 



H, - -hY,{S^ + S) + uY,iS] + SKS] + S-1) 
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(now assumed to be constituted by spins with S > 1/2), 
by proceeding along the same lines of appendix A. The 
resulting Hamiltonian Tls ■= (Cl^slO reads 

Hs = -h.Y.^^+^T. - l/25)(iD' + S/2] 

i i 

where h^, ^ h — U{2S — 1), and we have used the fact 
that, via a nontrivial calculation, one finds ((5^)^) — 
(1 - 1/25*) + S/2. Furthermore, stationarizing the 
action S = J dt{S\{ifidt — H^)]^) provides the equa- 
tions of motion for the variables L*, Li (see Eqs. (^), 
where LI is depending on L^, L* via the constraint intro- 
duced above for the spin expectation value components 
(if )2 -I- |LiP = 5^ Once the brackets 



(1 + 1?. 



\2\2 



2Sin 



dA dB 



dB dA 



(C7) 



have been defined, one can easily check that {L*,Lj} — 
2L'^/ih, {ij,i*} - L*/ih, and {i|,Lj} = -L^/ih con- 
sistently with ( p2[ ) , while the dynamical equations issued 
from the TDVP can be recovered as well from Ti.^ . Equiv- 
alently, the alternative form of the above Lie-Poisson 
brackets 



OA dB dB dA 



d^j dL'^ 



bj dLI 



(C8) 



can be reconstructed from (C7) when expressing i^'s as 



- {LIY e 



2 „'^<l>i 



(C9) 



through the action-angle variables 0i, if. This fact 
states at the classical level the equivalence between the 
HPR and the VR introduced in Sec. III. 



APPENDIX D: LOCAL PHASE DYNAMICS OF 
SHM HAMILTONIAN 

This appendix is devoted to calculate explicitly the 
form assumed by the spin Hamiltonian in proximity of 
the ground-state configuration in order to show how 
weakly excited states mimic the dynamics of the (clas- 
sic) phase model. Upon recalling that the ground-state 
configuration is characterized by 

Co = h*l{2U + qE,),^,^4>3. 
at each site j, considering the approximation 
BjL* + LiL* ~ 



2g cos{4>i — <pj) X 



where Pj = Lj — Co, g = — Co^ , and 



Co' 
2g 



^2 I I ^ 3 



2g 



leads to rewrite spin Hamiltonian (Q) as 

ns^c + Y,iup^+ uco' ~ KCo) 

i 

- gE, f{P,, P,) - gE,Y,cos{cj,, 

(iJ) ii,3) 



(Dl) 



Decoupling f{Pi,Pj) from the cosine term in the latter 
formula relies on the fact that {(j>j — (pi)'^ x f{Pi,Pj) is 
fourth order. The resulting model exhibits the QPM 
structure even if, within the present approximation 
scheme, the condition \Pj\ <C Co only concerns the spin 
dynamics. The geometry of the sphere (the spin con- 
figuration space) is involved instead when one imposes 
(3 = {"nil/S <^ 1 (m is the quantum number correspond- 
ing to Cq) in order to make explicit the (local) cylinder- 
like geometry characterizing the QPM in the equatorial 
region. 
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FIG. 1. The phase diagram of the BHM for D = 1. The error boxes are the QMC resu lts of Batrouni et al. in Ref. [9]. 
The dashed hnes are the result of the SCPE. Relatively to the first lobe {m = 1), Eq. ( |B24[ ) gives {tc/U) = 0.5. QMC gives 
{tc/U) = 0.43 ± 0.002 and the SCPE {tc/U) = 0.43. For {m = 3), QMC and SCPE give (tc/U) = 0.2 and (tc/U) = 0.18 
respectively. Our theory gives {tc/U) = 0.16. 



FIG. 2. The Phase Diagram of the BHM for D — 2 (continuous curve). The error box indicates the QMC tricritical point 
obtained by Krauth and Trivedi in Ref.[9]. For n = 1, {tc/U) = 0.25 while QMC gives {tc/U) = 0.244 ± 0.002 and SCPE 
provides {tc/U) = 0.272. 



FIG. 3. The phase diagram of the SHM for D — 1 and S = 10. Inside the lobes the ground state is paramagnetic. Elsewhere 
the system is in a canted phase. The lobes are obtained for m ranging from m = — S (corresponding to —2 < h/U < in the 
figure) to m = 1 (the upst lobe). The coordinates E^/U of the first two lobes' (corresponding to the —2 < h/U < 2 lobes) 
tricritical points decrease; they are "almost" constant for < h/U < 4 (named "inversion region" in the text). The SF phase 
is progressively reduced for increasing m {id. e. such that 4 < h/U < 16). 



FIG. 4. The phase diagram of the one dimensional S'-XXZ model for 5* = 10 and S = 15. The lobes are obtained for m 
ranging (—5* ... 1); in the figure, n = m + S ranges in {1 . . . 4}. We note that the SF phase is enlarged for increasing S. 



FIG. 5. The behaviour of the tricritical points {E^/U)tr as function of n for different values of S. We note that the region 
in which {Eg/U)tr is almost constant (named inversion region in the text) becomes wider for increasing 5*. Whithin this region 
the SHM's dynamics can be considerd as pure phase dynamics. This behaviour is named "QPM-like behaviour" in the text. 



FIG. 6. The phase diagram of the BHM compared with the phase diagram of the SHM for S — 55. The two phase diagrams 
coincide considering the zero-th order in a in formula (EU) with t — 2SE^ and h = /i 



FIG. 7. The phase diagrams of the SHM for 5* = 55 and /3 = 0.17 in D = 1. In Fig.6-a, the errorbars are the result of the 
QMC simulations of the QPM of Ref. [36]. 



FIG. 8. The phase diagrams of the SHM for 5* = 55 and /3 = 0.17 and for 5" = 100 and P = 0.22 in D = 2. The boxes are 
result of QMC simulation of Ref. [37] . 



15 



Fig.1 



5 - 



|i/U 3 - 



2 f: 



1 - 










1 


1 


1 


1 


1 1 












1D-BHM 






$ 






TDMFT _ 

SCPE 

- 


















^^^^^^^^^^^^ 















































1 


1 


1 




1 


1 1 







0.1 0.2 0.3 0.4 0.5 

t/U 



0.6 0.7 



Fig.2 



|i/U 




Fig3 




ID-SHM 
S=10 



_2 \^ I I I I I I I I I 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

E3/U 



Fig4 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 



0.02 
0.018 [- 
0.016 [- 
0.014 
0.012 

u 

^ 0.01 
0.008 
0.006 
0.004 |y 
0.002 




Fig.5 



20 



40 



60 



80 



r~ 

S=15 
S=20 
S=55 
S=100 



100 



n 



Fig.6 



|i/U 




Fig. 7 




